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Abstract, x decays into hadrons foresee the study of the hadronization of vector and axial-vector 
go : QCD currents, yielding relevant information on the dynamics of the resonances entering into the 

' processes. We analyse z — > 7t7t7tV T decays within the framework of the Resonance Chiral Theory, 

£Sj ! comparing this theoretical scheme with the experimental data, namely ALEPH spectral function and 

branching ratio. Hence we get values for the mass and on-shell width of the a\ (1260) resonance, 
and provide the structure functions that have been measured by OPAL and CLEO-II. 

> : 

o ■ 

cn INTRODUCTION 

o\ : 
o 

The implementation of Quantum Chromodynamics (QCD) in the energy region populated 
by light-flavoured resonances (M p < E < 2 GeV, being M p the mass of the p(770)) is a 
^3 ' demanding task that involves poorly known aspects such as bound and resonant states, 

duality and hadronization mechanisms. Though ad hoc Breit-Wigner parameterisations 
Qh| have been widely employed in the literature y],|2|] they are not necessarily consistent with 

the underlying theory, as they seem to violate the chiral symmetry of massless QCD 0,0]. 

T decays into hadrons allow to study the hadronization properties of vector and axial- 
vector QCD currents and, accordingly, to determine intrinsic properties of the dynamics 
generated by resonances [5]. At very low energies, typically E <C M p , Chiral Perturbation 
Theory (jPT) |@] is the Effective Field Theory of QCD. Still the x -» %%%v x decays, 
through their full energy spectrum, are driven by the p(770) and ai(1260) resonances 
mainly, in an energy region where the invariant hadron momentum approaches the masses 
of the resonances. Hence ^PT is no longer applicable to the study of the whole spectrum 
but only to the very low energy domain j2fl. The standard procedure that has been followed 
id to deal with these decays has been to modulate the amplitudes with a Breit-Wigner 
parameterization, fixing the normalization in order to match the leading 0(p 2 ) ^PT- 
Nevertheless, its deviation of the chiral behaviour at higher orders O,0] could spoil any 
outcome provided by the analysis of data. 

Lately several experiments have collected good quality data on x — > TtTZTtVx, such as 
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branching ratios and spectra or structure functions Their analysis within a model- 
independent framework is highly desirable if one wishes to collect information on the 
hadronization of the relevant QCD currents. 



THE RESONANCE CHIRAL THEORY OF QCD 

At energies E ~ M p the resonance mesons are active degrees of freedom that have 
to be properly included into the pertinent Lagrangian. The procedure, put forward in 
Refs. IKaTIiIi and known as Resonance Chiral Theory (R^T), is ruled by the approximate 
chiral symmetry of QCD, that drives the interaction of Goldstone bosons (the lightest 
octet of pseudoscalar mesons), and the SU (3)y assignments of the resonance multiplets. 
This construction is embedded within a comprehensive framework guided by the large 
number of colours (Nc) limit of QCD JH. The l/N c expansion tells us that, at leading 
order, we should only consider the tree level diagrams given by a local Lagrangian with 
infinite zero-width states in the spectrum. This is precisely the role of R^T. However in 
most processes, like hadron tau decays, we need to include finite widths that only appear 
at next-to-leading order in the large-Afc expansion and, moreover, we will only include 
one multiplet of vector and axial-vector resonances in our theory. Thus, in practice, we 
have to model this large-Afc expansion to some extent. 

The final hadron system in the % — > %%%v z decays spans a wide energy region, namely 
3/% < E < M x that is populated by many resonances. R^T is the appropriate framework 
to work with and we consider the Lagrangian ||4|,ll0(] : 

5 

+ -^kin + =^kin + 52 ^' ^VAP > (1) 
;=1 

where all the couplings are real, being F the decay constant of the pion in the chiral limit, 
and the operators & l VAP are given by : 

<p = (iV» v ,A, v ] X -) , 

^vap = i([V^ v ,A va ]h^} , (2) 

^vap = /([V% v ,A v > a > , 

<P = /<[V% v ,A a v K) , 

^vap = /<[V°%v,A' tv ]Ma> • 

The notation is that of Ref. ifioll . Notice that we are using the antisymmetric tensor 
formulation to describe the spin 1 resonances and, consequently, we do not consider the 
€?(p 4 ) #PT Lagrangian of Goldstone bosons liTTIl . 

Our Lagrangian theory has eight a priori unknown coupling constants, namely Fy, F A , 
Gy and A ; , i = 1,...5. Phenomenology could provide direct information on them; for 
instance Fy could be extracted from the measured r(p° — > e + e~), Gy from r(p° — > 



n + n~), Fa from T{a\ — > izy) and the Aj appear in r(ai — > p^r) or the T — > ;r7r7rv T 
processes themselves. It is conspicuous, though, that ^r^t is not QCD for arbitrary 
values of the couplings. Hence if we want to comprehend more about QCD in this non- 
perturbative regime we should try to learn about the determination of the couplings from 
the underlying theory II 1 311 . On this account we will implement several known features of 
the strong interaction theory in the following. 

The QCD ruled short-distance behaviour of the vector and axial- vector form factors in 
the large-Afc limit (approximated with only one octet of vector resonances) constrains the 
couplings of Jz^t in Eq. ([T]), which must satisfy ifllll : 

1-^ = 
F 2 

IFyGy - F y 2 = . (3) 

In addition, the first Weinberg sum rule, in the limit where only the lowest narrow 
resonances contribute to the vector and axial- vector spectral functions, leads to 

F^-Fl = F 2 . (4) 

In this way all three couplings Fy, Gy and Fa can be written in terms of the pion 
decay constant : Fy = VlF, Gy = F/y2 and Fa = F. These results are well satisfied 
phenomenologically and we have adopted them. In the next section we will comment on 
an analogous study of the A, couplings. 



THE AXIAL-VECTOR FORM FACTORS IN z~ -> {%%%)' v x 

The decay amplitudes for the t~ — > n + n^n^v r and t~ — > n~n°n°v T processes can be 
written as 

^± = - ^ V ud u v y ( 1 - y 5 K T ±ll , (5) 

T±n{phP2,Pz) = (MPi)MP2)n ± ( P 3)\^\0) 1 (6) 

as in the isospin limit there is no contribution of the vector current to these processes. In 
T±n(pi,P2,P3) the K + is the one in t~ — > iz + n~n~v z and n~ that in t~ — ► k~k°k°v t . 
The hadronic tensor can be written in terms of three form factors, F\, F2 and Fp, as lfl411 : 

= V?F\ + V?F 2 + Q tl F P , (7) 



where 



1 = V ~ ~W ) (pi_P3)v 



Q 2 

= p » +p » +p % . (g) 





(a) (b) 
FIGURE 1. Diagrams contributing to the hadronic amplitude 7±^ at &(p 2 ) %PT. 



The form factors F\ and F 2 have a transverse structure in the total hadron momenta 
and drive a J p = 1 + transition. Bose symmetry under interchange of the two identical 
pions in the final state demands that F\(Q 2 1 s 1 t) = F2(Q 2 ,t : s) where s = (p\+ p^) 2 and 
t = (P2 + P'i) 1 - Meanwhile Fp accounts for a J p = transition that carries pseudoscalar 
degrees of freedom and vanishes with the square of the pion mass. Its contribution to 
the spectral function of x — > nnnVt goes like rn^/Q 4 and, accordingly, it is very much 
suppressed with respect to that coming from Fi and F2. We will not consider it in the 
following. 

In the low Q 2 region, the matrix element in Eq. © can be calculated using #PT. At 
@{p 2 ) one has two contributions, arising from the diagrams in Fig. [2 The sum of both 
graphs yields 

T^ = T^{V 1 ,+V 2 ,} . (9) 

We now include the resonance-mediated contributions to the amplitude, to be evaluated 
through the interacting terms in J£r x t Eq. ©. The relevant diagrams to be taken into 
account are those shown in Fig. El We get 

T±^ = T ^f? V [a(£ 2 ,M)Vi M + a(£ 2 ,M)V 2M ] 

The functions a(Q 2 J s 1 t) and /3(<2 2 ,s,?) are : 

t ^r> s {IGy \ (2Q 2 — 2s — u u — s 1 

m ^ t) = -W + V)^ + F iQ ^J^ + FiQ \ t) -^ , (12) 



with 



2 
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F(£ 2 ,,) = -A -f + + (13) 
Q Or 




(a) (b) (c) 

FIGURE 2. Resonance-mediated diagrams contributing to T± jl . 



and depend on three combinations of the A; couplings in J^r^t* that we call Ao, A' and 
X" . Following the ideas outlined above we can get information on these combinations 
by implementing known aspects of asymptotic QCD. In particular we expect that the 
form factor of the axial-vector current into three pions should vanish at infinite transfer 
of momentum (Q 2 — ► oo). This is a consequence of the fact that its contribution to the 
spectral function of the axial- vector current correlator, being positive, has to add to other 
infinite hadronic positive contributions to reach the constant value evaluated within QCD 
iflfjl . Accordingly the proper behaviour of the T± il form factor imposes the constraints : 

2A' - 1 = , 

X" = . (14) 

Hence there is only one combination of couplings left unknown, namely Ao- 

Finally an additional comment on the result for 7±„ is required. The form factors in 
Eq. (TTOb include zero-width p(770) and ai(1260) propagator poles, leading to divergent 
phase-space integrals in the calculation of the x — > decay width as the kinematical 

variables go along the full energy spectrum. The result can be regularized through the 
inclusion of resonance widths, which means to go beyond the leading order in the 1 /Nc 
expansion, and implies the introduction of some additional theoretical inputs. This issue 
has been analysed in detail within the resonance chiral effective theory in Ref. [15] and, 
accordingly, we include off-shell widths for both resonances 0]. 



THEORY VERSUS EXPERIMENT 

To analyse the experimental data we will only consider the dominating J p = 1 + driven 
axial-vector form factors, that satisfy T +jl = —T-n hence providing the same predictions 
for both T~ — > Tl + Tl~TZ~ V z and T~ — > tz^tz q tz q v x processes in the isospin limit. 

We have fitted the experimental values for the T~ — > K + it^K^Vx branching ratio and 
normalized spectral function obtained by ALEPH [8] and we get a reasonable x 2 /d.o.f. = 
64.5/52. This is shown in Fig. |3j Hence we get the axial-vector ai(1260) parameters 
M A = (1.203 ±0.003) GeV and T ai (M 2 ) = (0.48 ±0.02) GeV, where the errors are only 
statistical. We also obtain a value for the still unknown combination of A ; couplings : 
Ao = 11.9 ±0.4. However, as pointed out in Ref. 11711 . this value seems too large when 
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FIGURE 3. Fit to the ALEPH data 0] for the normalized T~ -> % + %' V T . 



additional QCD constraints are imposed. The origin of the discrepancy could be the small 
sensitivity of the tau decay amplitude to this parameter, as it only appears multiplied by 
the mass of the pion (fT3"l) . together with an improvable implementation of the off-shell 
width of the ai(1260) resonance. 

Ultimately we predict the integrated structure functions wa, wc, Wd and we [ 14], that 
we compare with the experimental results for t~ — > n~n n°v T in Fig.|4] In spite of the 
large errors the predictions follow notably the depicted trend. As analysed in Ref. 0] a 
variance between the ALEPH data on one side and the OPAL and CLEO-II on the other 
are at the origin of the seemingly inconsistent result for wa in the high Q 2 region. 

In conclusion it can be inferred that, within the present experimental errors and for the 
studied observables, there is no evidence of relevant contributions in x — > %%%v x decays 
beyond those of the p(770) and ai (1260) resonances. 



ACKNOWLEDGMENTS 

J.P wishes to thank Pietro Colangelo and Fulvia De Fazio for the remarkable organiza- 
tion of the QCD@Work 2005 meeting in Conversano (Italy). This work has been sup- 
ported in part by MEC (Spain) under grant FPA2004-00996, by Generalitat Valenciana 
(Grants GRUPOS03/013, GV04B-594 and GV05/015) and by HPRN-CT2002-003 1 1 
(EURIDICE). 



2e+05 



a 



4x10 



2x10 




Q* (GeV ) 



Q" (GeV ) 



le+05 




2e+04 



-5e+04 




-2e+04 



Q 2 (GeV 2 ) 



Q 2 (GeV 2 ) 



FIGURE 4. Theoretical values for the wa, wc, wd and we integrated structure functions in comparison 
with the experimental data from CLEO-II (solid) and OPAL (dashed) [9]. 
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